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ABSTRACT 

We propose to use a simple observable, the fractional area of "hot spots" in weak gravitational lensing mass 
maps which are detected with high significance, to determine background cosmological parameters. Because 
these high-convergence regions are directly related to the physical nonlinear structures of the universe, they 
derive cosmological information mainly from the nonlinear regime of density fluctuations. We show that in 
combination with future cosmic microwave background anisotropy measurements, this method can place con- 
straints on cosmological parameters that are comparable to those from the redshift distribution of galaxy cluster 
abundances. The main advantage of the statistic proposed in this paper is that projection effects, normally the 
main source of uncertainty when determining the presence and the mass of a galaxy cluster, here serve as a 
source of information. 
Subject headings: 



1. INTRODUCTION 

Weak gravitational lensing (WL), i.e., the coherent 
distortion of images of faint distant galaxies by the 
gravitation al tidal field of the i nterve ning matter dis- 
tribution dTyson. Wenk. & Valdesl 1 1990h . has been es- 
tablished as a powerful co smological tool [see, e.g., 
iBartelmann & Schneider! d2001l) for a review]. Since this ef- 
fect is purely gravitational, it directly probes the matter dis- 
tribution along the line of sight, thus providing a way to un- 
derstand the nature and the evolutionary history of the uni- 
verse that is relatively insensitive to how light or baryons 
trace dark matter. Since the earlie s t measurements of WL 
by ga laxy clusters dFort et al.l 119881 : iTyson. Wenk. & Valdesl 
1990), WL by large-scale structure, know n as "cosmic 
shear ", has been detected based on optical dWittman et alJ 
2000; IVan Waerbeke et all l200(A iBacon. Refregier. & ElTisI 
2000; Kaiser, Wilson, & Luppino 2000) and radio observa- 
tions dChang. Refregier. & He lfand 20()3)- 

Future WL surveys, such as the Large Synoptic Survey Tele- 
scope (LSST) 5 , will measure the cosmic shear field with great 
precision over half the sky. Recent attention has focused on 
how best to extract information from the shear/convergence 
field in order to constrain cosmological parameters. The pri- 
mary focus so far has been us i ng the standard two-poi n t statis- 
tics dBlandford et all Il991h iMiralda-Escudl I199U iKaiserl 
119921) . which probes the underlying matter power spectrum 
in projection. In this paper, we propose to use a simple, direct 
observable: the fraction of high signal-to-noise ratio (S/N) 
points detected in WL surveys, as another discriminator of 
cosmology. We utilize the results of an N-body simulation to 
quantify both the theoretical predictions and the observational 
uncertainties of this statistic. 

As is well-known, the common two-point statistics do not 
contain all the statistical information of the WL convergence 
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field, as the nonlinear gravitational instability induces non- 
Gaussian signatures in the mass distribution and hence in the 
WL convergence field as well. Unfortunately, there is no com- 
plete statistical analysis in practice for the underlying mat- 
ter density field in the nonlinear regime. Previous studies of 
weak lensing statistics on small angular scales, correspond- 
ingly, were restricted to the low-orde r statis tics, using the 
halo model [see, e.g.. ICoorav & Shethl d2002l) for a review], 
or the "scaling ansatz" dHamilton e t al. 19911 ), later extended 
and calibrated by N-body simulations dJain. Mo. & White! 
11995b iPeacock & Doddsl 119961: ISmith et alJ 120031) . for the 
nonlinear evolution of clustering. These low-order statis- 
tics include: the two-point correlation funct ion or equiv- 
alentl y the power spectrum, at small scales (I Jain & Seljakl 
1 1997b , the three-point correla t ion fu nction or the bispec- 
trum dTakada & Jain] l2003albl 120041). the third-order mo- 
ment: skewne ss (Bernardea u. van Waerbeke. & Mellieril997j 
Uain & Seljakl 119971: iHuil Il999f) and the fourth-order mo- 
ment: kurtosis (Takada & Jainl [2002b . An alternative ap- 
proach is to use the redshift distribution of nonlinear object 
abundance (Hai man. Mohr. & Holder] 1200 ll). such as shear- 
select e d galaxy clu s ter samples ([W einberg & Kamionkowski 
120031: IWang et alJ l2004t iFang & Haimanl |2007l). using 
the Press-Schech ter prescription dPress & Schechterl 1 19741 : 
iBond et aTJ Il991h with calibration by N-b ody simulations 
dSheth & Tormenll 19991: 1 Jenkins et al]|2001l). In addition, the 
"ratio-statistic" ([Jain & Ta ylor 2003; Bernstein & Jainl 120041: 
[ Song & Knoxl2004l:lHu & Jainl2004t[Zhang. Hui. & Stebbinsl 
2005) has been constructed to use the geometrical informa- 
tion from tomographic shear/convergence power spectra on 
small scales. These analyses have shown that comparable in- 
formation may be contained in the linear and in the nonlinear 
regime. 

The statistics mentioned above are by no means a full 
characterization of the convergence field, yet they all face 
challenges, either from the theoretical or the observational 
side, which need to be addressed. Including higher-order 
statistics might give a substantial increase in information, 
but they are in practice noisy and computationally inten- 
sive. To take advantage of the synergy between these statis- 
tics, one needs to account for their covariance, which is dif- 
ficult to calculate analytically. One hope is to run large 
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simulations and directly measure this covariance. Galaxy 
cluster samples selected optically, by their X-ray flux or by 
their Sunyaev-Zel'dovich effect signatures, are limited by the 
uncertain astrophysics when modeling t he mass-observable 
relations and require "self-calibration" (jMa iumdar & Mohr 
l2003HWang et al.ll2004t Lima & Hull2005l) . as the mass func- 
tion is exponentially sensitive to errors of limiting mass. 
Shear-selected galaxy clusters, on the other hand, have the 
advantage that the selection function can be determined 
ab initio by N-body simulations, since only gravity is in- 
volved. However, proje ction effects result in false de- 
tections, missing clusters (IWhite. v an Waerbeke, & Mackevl 
2002|:lHamana. Takada. & Yosh ida 2004; Hennaw f&Spergell 



2005) , and producing significant uncer tainty in the cluster 
mass derived from the shear signals (M etzler et alj 119991: 
Ide Putter & Whi te 2005), degrading the cosmological infor- 
mation content. 

In this paper, we instead focus on the one-point probabil- 
ity distribution function (PDF) of the WL convergence field. 
There are three main motivations. First, the one-point PDF 
is a simple yet powerful tool to probe non-Gaussian features, 
and since the non-Gaussianity in the convergence field is in- 
duced by t he growth of structure, it holds cosmological in- 
formation ( Reblinsk y et alJl999tlJain. Seliak. & Whitell2000t 
iKruse & Schneiderl2000tlValageas. Munshi. & Barberl2005l) . 
These previous works have shown that the one-point PDF 
is capable of discriminating cosmologies with different fi m , 
such as an open cold dark matter (CDM) model, a flat cosmo- 
logical constant-dominated (A)-CDM model, and a standard 
CDM model. 

Second, and the primary motivation of this work, is that the 
fractional area statistic we propose in this paper takes into ac- 
count projection effects by construction. Determined by the 
high-convergence tail of the PDF, it is an analog of Press- 
Schechter formalism, thus similar to the abundance of galaxy 
clusters but without contamination due to projection effects. 
This statistic also utilizes information mainly from the nonlin- 
ear regime and complements the well-established statistics in 
the linear regime. The goal of our work is to give a more 
quantitative assessment of the statistical information in the 
nonlinear regime (particularly focusing on the properties of 
dark energy) provided by this fractional area statistic and its 
complementarity to other probes of cosmology, e.g., the cos- 
mic microwave background (CMB) anisotropics. 

Third, besides utilizing the cosmic shear field 
(Zhan g~& Pen! [2005b which we will focus on in this pa- 
per, there are several other observational techniques which 
could be used to map out the convergence PDF. For example, 
with forthcoming large samples of high redshift supernovae 
from LSST-like or Joint Dark Energy Mission (JDEM)-like 
surveys, one could measure the magnification distribution 
of these standard candles due to lensing by the large scale 
struct ure in the foreground, and construct the convergence 
PDF (iDodelson & Vallinottoll2006t ICooray. Holz & Hutererl 

2006) ; it is also possible to measure the conve rgence fiel d 
through the statistics of cosmic magnificatio n dJainl 1 2002), 
using, for instance, 21 cm-emitting galaxies (IZhang & Pen! 

20051 EqqS)- 

The rest of this paper is organized as follows. Our ba- 
sic calculational methodology, including a calibration of the 
(co)variance using simulation outputs, is described in §[2] Re- 
sults of using the fractional area statistic for an LSST-like WL 
survey are presented in § [3] with discussions of its comple- 
mentarity to other dark energy probes, as well as of various 



uncertainties. Conclusions and implications of this work are 
given in § |U Finally, in a series of appendices, we show the 
details of our calculations. 

2. CALCULATIONAL METHOD 

2.1. Convergence Field with Gaussian Smoothing 

Consider some source galaxies detected in a WL sur- 
vey, with a redshift distribution of n(z) and a surface 
density of n g = J n(z)dz- The shear field is mea- 
sured from the distortion of their images. A conver- 
gence field, or sometimes called mass map, can be 
recon st ructed from the shea r field ((kaiser & Squires! 
1993t [Squires & Kaiserl [19961 iBartelmann et all 119961: 
Van Waerbeke. Bernardeau & M ellier 119991) . smoothed 
over scale Qq with a Gaussian window function 
W G (9) = exp(-e 2 /0 2 c )/(ir0 2 c ). This map is a sum of the 
true, smoothed convergence field and the noise field due to 
the randomly oriented intrinsic ellipticities of source galaxies: 



K(6) = K (S) (e) + K(N) (6). 



(1) 



Both fields are assumed to have zero mean and are statistically 
isotropic (ensemble average is same along each line of sight). 

Under the assumption that the correlation of intrinsic el- 
lipticities and the clustering of the source galaxies can 
be neglected, and n o other systematic errors are present, 
IVan Waerbekd (|2000) has shown that the noise field r^) can 
be modeled as a Gaussian random field with variance: 



4tt6 



(2) 



where a t is the root-mean-square value of the intrinsic ellip- 
ticity of source galaxies. The two-point correlation function 
of K(N) induced by the smoothing is 



C (N 0) = a (N) exp[-— T ) , 



(3) 



where 6 = \6\ —6z\ is the angular separation between two 
points in the field. 

The true convergence field, on the other hand, is essentially 
non-Gaussian. Its one-point PDF is skew, because K($) has a 
minimum value, corresponding to an totally empty path be- 
tween the source and the observer. The convergence without 
smoothing, by using Born approximation, is a weighted pro- 
jection along a particular line of sight of the mass density pe r- 



turbation field (Bernardeau, van Waerbeke, & Mellier||i997|) 



k(9)= / d X WUx)5(x9,x), 



(4) 



where \ is the comoving distance to redshift z, and S is the 
over-density at comoving distance \- Taking 5 = -1 every- 
where along the line of sight, gives the minimum value of 

K(S)- 



K(S)mm — »vmin ' 



<^XWlen(xX 



(5) 



as the smoothing has no effect on a constant 6 . For a fixed 
source plane at redshift z s corresponding to n(z) = n g 5(z-z s ), 

6 Strictly speaking, this is not true if one measures the minimum value 
from simulations or real surveys, due to the finite-sampling effect: it is more 
rare to have totally empty lines of sight for the whole smoothing aperture. 
This effect thus results a higher "K m i„" for a smoothed field. 
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the weight function is given by (assuming a fiat universe for 
simplicity): 

wux,xs)=ln m (jf) a+z)x(i-^-)n(x s -x), (6) 

where Xs = Xfc) denotes the comoving distance to the source 
plane and TC is the Heaviside step function. Note that this 
quantity k^, which is a source of the non-Gaussianity of the 
convergence field, depends on the source redshift and on cos- 
mology. The covariance and the correlation function of the 
smoothed convergence field, using Limber's approximation, 
are given in Appendices A and C as Eq. (ICU and iA3i . 

Points with high S/N in the convergence field are of both 
experimental and theoretical interest: not only are they related 
to the underlying physical density field at high fidelity, but 
they capture the important features of the nonlinear regime 
as well. Mathematically, they correspond to the so-called 
"excursion set", whose properties have been studied exten- 
sively in other similar cont exts (most famou sly to describe 
the halo mass function, e.g. jBond et al.llT99 lb . For pedagog- 
ical purposes, let us consider the excursion set E v of the true 
smoothed convergence field at redshift z s , defined as the union 
of all points with K(s) > ^cr(N) (see §[3]for including the noise 
field). The fractional area of E v is then given by 



1 



F(v,z s ) = — = — Id 9H[k(s)(9,z s )-v<J(n)], 



(7) 



where A v is the area of the excursion set E v , and A tot is the 
total sky coverage of a survey. K( S ) is a statistical quantity, so 
the fractional area has a mean value of 



dK(S)Pl(K(S)), 



(8) 



where Pi is the one-point PDF of K(S). The spatial integra- 
tion cancels with A tot because of the statistical isotropy of 
the field 7 . Similarly, the fractional area has a variance de- 
termined by the two-point PDF, which is given in Appendix 
CbyEq. CD. 

2.2. Statistical Properties of the True Convergence Field 

Let us first focus on the one-point PDF of K(s) at dif- 
ferent redshifts: Pi(k<js),Z s ), which determines the mean 
fractional area {T{v,z s )). Several wo rks have used the 
"stable-clustering ansatz" (Peeblesl 119801) t o derive the one- 
point PDF of the true convergence, with dValageasI [2000b; 
Muns hi & Jainl2000t IWang. Holz. & MunshifcOOl and with- 
out ( Valageas 2000a) smoothing, calibrated by N-body simu- 
lations. The important conclusion of these works is that there 
exists a universal one-point PDF, well approximated by a two- 
parameter family: the variance of the reduced convergence 
( K ts)) I K min an d me minimum convergence K m i n . 

In this section, we concentrate on the constraint equations 
of the one-point P DF required by the no rmalization, the mean, 
and the variance dDas & O striker 2006). By normalizing K(S) 



to K 



: K (S)/ 



, the constraints can be written as 



cIk'Pi(k') = 1, 



7 Equivalently, we can assume ergodicity, so that the spatial average is 
equal to the ensemble average. 

8 These three constraints will be sufficient to specify a functional form 
of PDF with no more than three fitting parameters. For more complicated 
models, one needs to go beyond the variance and specify for the higher-order 
moments their dependence on the variance and the minimum. 



/oo 
dK'Pi(K')K' = 0, 



(9) 



We consider three different fitting formulae for the one-point 
PDF of the convergence field from the literature. In the next 
section, we will use WL simulation outputs to assess the ac- 
curacy of the high convergence tail of these three models. 

The firs t model we conside r is the log-normal distribution 
studied in Taru ya etail d2002), and is given by 



Pl(K(S))dK(S) = 



dn ( s) 



(10) 



V27 



: exp 



ZjV 



[ln(l + K W /|/t m in|) + gL/2] ; 
2(t In 



Only one parameter, a]^ = ln(l + (rtm) / K 2 ^), is needed to 
satisfy the constraint equations. 

The seco nd model is the stretched Ga ussian distribution 
proposed in lWang. Holz. & Munshjl (120021) . 



dK(S) 

Pi(K(S))dK(s) = | rC norm exp 



(l+K( S )/|K min |)-77p eak 



p(l + K(,S)/|Knija|)« 

(11)' 

Here we closely follow the original notation, except the 
PDF is now written in terms of n^s) instead of rj = 1 + 
K(S)/\Kmin\- As stated earlier, the four fitting parameters C n0 rm, 
>7 P eak, P an d q depend only upo n (k^) / [see Eq. (8) 
of IWang. Holz. & Munshil (l2002l) l. The complication of this 
model is that there is a need to impose an upper limit of the 
convergence K max in order to get the correct variance. 

The third model takes t he form of a modified log-normal 
distribution as proposed in lDas & Ostriker! ([2006): 



Pi(KQS))dK($ = 
xNexp 



dn (S) 



, K m in | + K(S) 

[ln( 1 + K min I ) + S 2 /2] 2 [ 1 + A/{\ + K (S) / 1 K^n |)] 

2E 2 

Given (k 2 s) ) /k^, we can uniquely specify N, £ and A by 
numerically solving the three constraint equations. The de- 
pendence of t hese three parameters on the variance is plotted 
in Figure 2 of iDas & Ostriker! d2006l) . 

For the analysis in this paper, we emphasize that the cosmo- 
logical dependence of these PDFs, i.e., information on dark 
energy, enters only through the variance of the reduced con- 
vergence 9 : (rfs)) / K min' an( ^ through the minimum: « m ; n as an 
additional scaling factor. 

Once the form of the convergence PDF is specified, the 
mean fractional area is determined by its integral [Eq. {§)]. 
The (co)variance of the fractional area for various thresholds 
and redshifts, Cov[J r (iJ,,Zi),^F(i',Zj)], on the other hand, re- 
quires knowing the joint two-point PDF of the convergence 
field, as we can see, e.g., in Eq. CD- To the best of our knowl- 
edge, among the three cases considered above, the two-point 
function can be computed analytically only in the log-normal 

9 The universal halo mass function similarly depends on cosmology only 
through cr(M,z), the rms value of the matter density fluctuations. However, 
one should keep in mind the essential difference that here the convergence 
variance is calculated using the nonlinear power spectrum, whereas for the 
cluster mass function, the variance is calculated using the linear matter power 
spectrum. 



(12) 
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model; the results are presented in Appendix C. Note that we 
do not use the covariance matrix of the log-normal model for 
our actual calculations - we will instead utilize simulation 
outputs to directly measure the (co)variance. The results for 
the log-normal model will be used only for comparison, and 
to justify extrapolations from simulations, as we discuss in the 
next section. 

2.3. Simulation and Model-Fitting 

The simulation outputs we use are those of lWhitel d2005l) . 
There are a total of 32 convergence fields from two indepen- 
dent N-body simulation runs. Each field is 3 x 3 deg 2 , di- 
vided into 1024 x 1024 pixels, so that the angular size of a 
pixel is ^10.5 arcsec. There are three source redshift planes 
at z s = 0.6, 1.1, and 1.9 for each field, representing the mean 
redshift of three source galaxy bins of an LSST-like sur- 
vey. The cosmological parameters are adopted from first- 
year measurements by the Wilkinson Microwa ve Anisotropy 
Probe (WMAP) 10 , as summarized in Table 1 of ISpergel et al.l 
(2003): a spatially flat ACDM model with a scale-invariant 
initial scalar power spectrum (n s =1) and present-day nor- 
malization erg = 0.9. The matter density is Vt m = 0.28 and the 
baryon density is il/, = 0.049. The Hubble constant is Hq = 70 
km/s/Mpc. 

We first apply a Gaussian window function to smooth all 
the fields, with 6c chosen to be 1 arcmin (~ 5.7 pixel): 



K(.S)(i,j)= E K(i+m,j + n) 



exp 



n 2 +n 2 <32 2 



5.7 2 



E eX P 

m 2 +n 2 <32 2 



5.7 2 



(13) 



(i, j) and (m,n) are integers used to label the coordinates. For 
each pixel, we smooth it with nearby pixels within a circle of 
a 32-pixel radius which is roughly 6 arcmin, as the Gaussian 
weighting is negligible outside. We discard 32 pixels on each 
side of the field to avoid possible edge effects, and use only 
960 x 960 pixel 2 in the center of each smoothed field for the 
analysis. 

In Figure Q] we show fits to the normalized histograms of 
the smoothed convergence fields for the three different PDFs 
discuss above (using a bin-size of A/t(j) = 0.0025). As de- 
scribed earlier, for all three models, there are only two pa- 
rameters to fit: the minimum of the distribution and the vari- 
ance. For the log-nor mal model and the model proposed by 
IDas & Ostrikerl (120061) . all fitting parameters are fixed once 
these two quantities are given. H o wever , for the model pro- 
posed by Wan g. Holz. & Munshi (2002), there is a further 
complication in determining K max . For this model, we have 
therefore relaxed the assumption above, and allowed all pa- 
rameters to vary simultaneously, except we imposed a prior 
that they lie in the ranges r oughly matching Figure 2 in 
IWang. Holz. & Munshil(l2002l) : 0.6 < r; peak < 1, < w < 0.5, 
and 0.5 < q < 2.5. 

Recall that the log-normal distribution would be a sym- 
metric parabola on a log-log plot. The actual distribution 
from simula tions, however, is skew ed toward the high-K(s) 
tail [see also lDas & Ostrikerl d2006l) l. so the log-normal PDF 
under-predicts the tail of the distribution. In contrast, the 
Wang. Holz. & Munshil (120021) model over-predicts the distri- 
bution in this range, as the distribution is cut at K max . Overall, 

10 map.gsfc.nasa.gov 



the model bv lDas & Ostrikerl ([2006 ) goes through the simula- 
tion data points remarkably well at all three redshifts 11 . 

A sl ightly different expression for K m ; n was given in lLinderi 
d2007l) . taking into account density inhomo geneities by ap- 
plying an e xtended Dyer-Roder formalism (Dyer & Roeder 
19721 Til 9731) to a compl etely empty line of sight [also see 
Seitz & Schneider ( 1992)]. The correction to the commonly 
used Eq. © is only large when the sources are at very high 
redshifts (10% for our fiducial cosmology at redshift z = 3). In 
principle, K m ; n affects the high-tail of the distribution, since it 
is used to normalize the universal PDF. One could also di- 
rectly measure this minimum quantity from the simulations. 
However, the low end of the PDF might be really suppressed, 
because a completely empty line of sight is extremely unlikely 
in a finite-size simulation; the smoothing procedure makes it 
even less probable. We have found that the minimum conver- 
gence measured from simulations is neverthe less larger tha n 
theoretical predictions, either using Eq. (| 5j) or [Linden d2007l) . 
The best-fit K m ; n of the IDas & Ostrikerl diOQq) model, how- 
ever, agrees well with the measured minimum value. This is 
a ~ 20-50% effect, depe nding on the source redshift [see a 
more detailed discussion in lTaruva et al.l d2002l) l . One intrigu- 
ing question is which, among the minimum values mentioned 
above, could lead to a more universal PDF? Particular, for the 
purposes of this paper, which one leads to the most univer- 
sal high-tail distribution? We hope to investigate this question 
in future work, using larger-volume simulations with various 
cosmologies. In this paper, we simply use Eq. (0) to calcu- 
late the derivatives of the fractional areas with respect to the 
cosmological parameters. 

To conclude, for the purpose of taking derivatives in the 
Fisher matrix analysis (see below), we adopt the universal 
one-point PDF provided by lDas & Ostrikerl (120061) . with the 
variance of the convergence field co mputed by using the 
nonlinear matter power spectrum from Smit h et alj d2003l) as 
given by Eq. dCU and the minimum as given by Eq. (0). 

Next, we measure directly from the simulation outputs the 
covariance matrix of the fractional areas with various thresh- 
olds and redshifts: Cov[J 7 (fj,,Zi),J r (i',Zj)]- The variance is 
given by 



Var[^,z s )]: 



N f -1 



E 



N„ 




f -,N( m Hv,z s ) 



N„ 



.(14) 



Here M m) (^,z,) is the total number of pixels with signal above 
the threshold v<J{n), in the m-th field at redshift z s \ Nf = 32 is 
the total number of fields and N p = 960 2 is the total number of 
pixels in each field. Similarly, the covariance is given by 



Nf-l 



(15) 



N f 

E 

i n= i 



N (m) (^z i )N (m Hi',Zj) 1 



N„ 



N„ 



3 m=l 



m ) 



1 1 One might wonder how one curve in Figure \T\ can lie entirely below 
another, when the total area under each of the curves is normalized to unity. 
Notice that the logarithmic y-axis shows a range of many order of magnitudes 
- the tails, where the three models differ, visually dominates these figures, 
but the area under these tails is negligible compared to the total area of the 
distribution, which lies around the slightly offset peak of each PDFs. 
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To extrapolate the results we obtain from a 9 deg 2 simula- 
tion field to a 20,000 deg 2 LSST-like survey, we also measure 
the covariance matrix as a function of the field size. To do 
this, we divide each field into several smaller sub-fields. The 
following cases are considered: N p = 480 2 , 240 2 , 120 2 , 60 2 , 
and 30 2 . For each case, we fix Nf = 32. For example, when 
N p = 480 2 , we first use the upper left 480 2 pixel sub-field of 
all 32 fields to measure the covariance using Eq. (116b . We 
repeat the same procedure for the other three sub-fields, and 
then take the average of the four results. 

In Figure |2] we plot the (co)variance of the fractional area, 
as a function of scale, for different thresholds and redshift bins 
measured from the simulation outputs. For small fields, we 
effectively have more realizations to measure the covariance 
matrix, for which we plot the average. For comparison, we 
also plot, as dotted curves, the values calculated using the log- 
normal model. 

One expects that when the size of the field mdx becomes 
larger than the correlation length of these "biased regions" in 
the convergence field, the (co)variance scales with respect to 
the survey size at least as steeply as oc A^ t (coinciding with a 
Poisson-like scaling; see Appendix C). On smaller scales, the 
slope is generally shallower due to correlations. For example, 
we can observe this trend analytically in the log-normal cas e 
- the covariance is given by a Poisson expression [Eq. ( IC9H 
for large scales and an additional non-Poisson correction term 
[Eq. (IC11H , The simulation points in Figure [2] show that for 



low-redshift bins, the scale of 1024 pixel (3 deg) is already 
within the Poisson region. For higher redshift bins, the slope 
at 1024 pixel is shallower than Poisson. However, the slope 
does agree well with the log-normal model, which predicts 
that the scale of 1024 pixel is very close to the edge of the 
Poisson region. Thus, for our Fisher matrix analysis, we esti- 
mate the covariance matrix of the fractional area for an LSST- 
like survey by Poisson-scaling the measured value from the 
simulation outputs with the largest field size of N p = 960 2 
pixel. This extrapolation to 20,000 deg 2 is shown explicitly 
in Figure |2] as the solid curves. The justification for this ex- 
trapolation, using the log-normal model, is admittedly still 
somewhat heuristic, but measurements of the covariance out 
to larger scales require larger size simulations; we defer this 
to future work. 

2.4. Random Ellipticity Noise 

In order to incorporate the presence of noise into our fore- 
casts, we start with the simplest case: modeling the noise due 
to the intrinsic ellipticity as a two-dimensiona l Gaussian ran- 
dom field that is uncorrected with the signal (IVan Wae rbeke 
2000; |Mn_&vanWaerbeke 2000). We use a Gaussian ran- 
dom number generator with zero-mean and a variance of (j(M) 
to generate noise fields and co-add onto the original conver- 
gence fields. We adopt a e = 0.3, 9c = 1 arcmin, and n g = 40/3 
arcmin" 2 (the denominator "3" reflects the redshift binning) in 
Eq. (O for an LSST-like survey, thus cr^j ~ 0.023. We then 
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FIG. 2. — The top left panel shows the variance Var[J-"(/x = 2.z,)] within three distinct redshift bins; the top right panel shows the covariance Cov[jF(^i = 
2,Zi)-,T(v = 2.5, z,)] within the same three redshift bins. The bottom panel shows Co\[^F(/i = 2,zi),T{v = 2,z.j)] across three different pairs of redshift bins. All 
symbols in each panel show results measured directly from the simulation outputs. The solid curves show the estimate of the covariance matrix through Poisson- 
extrapolation, which we use in the actual calculations. The dotted curves are calculated using the log-normal model, which are plotted only for comparison. 



repeat the same Gaussian smoothing procedure and measure 
the covariance matrix of the fractional area. In this case, the 
one-point PDF is simply a convolution of the original distri- 
bution and the Gaussian PDF. 

Modeling of the noise field as above might be overly sim- 
plified, even in the c ase when all other sys tematics are ab- 
sent. This is because IVan Waerbekd (|2Q00) only worked in 
the WL limit, in which the estimator of the convergence is 
safely linear. We here typically use a signal threshold of sev- 
eral times <T(jV), reaching values of order ~ 5- 10%. Con- 
struction of the convergence field from the observed reduced 
shear is nonlinear, and therefore gets non-negligible higher- 
order corrections, which will complicate the expression for 
the convergence error. Besides the intrinsic ellipticity, there 
will be other systematic effects, due to the point spread func- 
tion, mirror distortion, etc, which will increase the scatter, and 
may produce an unknown bias in measuring the cosmic shear. 

In our analysis, we adopt an approach similar to what has 
been done in the literature for cluster counts in order to ac- 
count for the uncertainties of the mass-observ a ble relations 
(Maiu mdar&Mohrj 120031: iWang et alJ 12001 iLima & Hul 
120051) . We still make the simple assumption that the over- 
all convergence noise is Gaussian and uncorrected with the 
signal, but with unknown variance and bias, which will be in- 
cluded as additional nuisance parameters that are fit by the 
survey itself, simultaneously with the cosmology parameters 



("self-calibration")- We assign the variance and bias in each 
redshift bin as free parameters, and their fiducial values are 
{(j(Ai), Kbias} = {0.023,0}, chosen to be same for all redshift 
bins. For WL powe r spectrum and bispectrum tomography, 
iHuterer et al.1 (|2006) have shown that there is sufficient in- 
formation in these statistics themselves to successfully self- 
calibrate additive and multiplicative shear systematics. As we 
will see below, the nuisance parameters we consider here like- 
wise do not significantly degrade the cosmological constraints 
derived from the one-point statistic, provided we have some 
external knowledge of what the intrinsic ellipticity distribu- 
tion is. The reason for this is essentially that the PDF contains 
both the tomographic and the shape information (similarly to 
the case of cluster counts, where the nuisance parameters can- 
not simultaneously mimic changes in both the shape and the z- 
dependence of the mass function caused by variations in cos- 
mology). 

2.5. Error Forecasts 

The usual Fishe r matrix technique 

(Tegmark. Tayl or? & Heavens! Il997t) is employed to as- 
sess the cosmological sensitivity of the fractional area 
statistic. This method allows a quick exploration of the 
parameter space and gives a lower bound to the statistical 
uncertainty of each model parameter to be fit by future 
experimental data. 
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We consider a WL survey with specifications similar to 
that planned for LSST: a sky coverage of A tot = 20, 000 deg 2 , 
source gal axies with the redsh ift distribution adopted from 
Eq. (18) in Son g & Knoxld2004l) . normalized to a total surface 
number density of n g = 40 arcmin" 2 , and the intrinsic elliptic - 
ity dispersion a t = 0.3. The galaxies are divided into three 
bins with mean redshift z s = 0.6, 1.1 and 1.9 so that each bin 
contains the same number of galaxies. The Gaussian smooth- 
ing is taken over the scale Qq = 1 arcmin within the nonlinear 
regime. Seven different S/N thresholds, from v = 2.0 to 5.0 
in increment of Av = 0.5, are considered simultaneously for 
each redshift bin to utilize the information contained in the 
shape of the PDF (see §[5] below for discussions of different 
choices of smoothing scale and S/N thresholds). 

The full Fisher matrix is g iven by Eq. (15) of 
Teg mark. Taylor. & Heavens! (119971) : 



fa/3 



= I T r 

2 



C C Q C l Cj3 + C 1 [T^aT'p + TjjT', 



(16) 

where T is the list of mean fractional areas above different 
thresholds and at different redshifts, arranged into a vector, 
and C is the covariance matrix between them. We have used 
the standard comma notation for derivatives. 

In this paper, we neglect the cosmological dependence of 
the covariance matrix, i.e. the sample variance 12 . It is 
straightforward to spell out the last two terms, which gives 



4t,[c 



dp a 



(17) 



where p a includes the cosmological parameters of interest, 
as well as possible nuisance parameters. The set of different 
redshift bins and S/N thresholds are labeled by i or j, and /i 
or v, respectively. 

The covariance matrix C is estimated by extrapolation from 
simulations, as discussed above. Note that in Eq. ( TP7I ). this 
covariance matrix is needed only for the fiducial cosmology. 
The subscript "(jj,,i)(v,j)" denotes the corresponding matrix 
element. The choice of cumulative S/N thresholds instead 
of independent bins is made to simplify the calculations in 
the Appendices. Note that the off-diagonal terms of the co- 
variance matrix, reflecting the cross-correlations between dif- 
ferent redshift bins and S/N thresholds, are nonzero so that 
double-counting information from overlapping lens redshift 
and S/N ranges is avoided. 

We consider a spatially flat cosmological model with seven 
parameters. The fiducial values are close to those adopted 
in the simulations: {iloE,wo,w a ,fl m h 2 ,ilbh 2 ,(Js,n s } = 
{0.72,-1,0,0.137,0.024,0.9,1} (note a higher cr 8 than 
WMAP three-year result). The dark energy equation of state 
is parametrized as dChevallier & Po larski 200 It lLinderll2003l) 

z 

w(a) = wo + w a (l-a) = wo + w fl - . (18) 

l+z 

CMB anisotropies have played a crucial role in firmly 
establishing the current "standard model" of cosmology 

12 Just like the full Fisher matrix, the sample variance matrix is positive 
definite, because it is a special case when the observables have mean val- 
ues of, for example, zero. Thus we expect better constraints if this term 
was included, since it contains information from two-point statistics of these 
WL "biased" regio ns. This is again analogous to the case of galaxy clus- 
ters, where Lima & Hu 12005) have shown that sample-variance helps with 
self-calibration. We hope to explore this in a future work. 



dSpergel et all I2003L l2007h . To take advantage of the fu- 
ture high precision data which will be available from planned 
CMB anisotropy measurements, a survey with specifications 
similar to the Planck surveyor (Planck) 13 is considered. The 
Fisher matrix for the CMB temper ature and polar ization 
anisotropies is constructed as given in Wa ng et ail (l2004f) . 

3. RESULTS AND DISCUSSION 

Before we present the results, it is useful to keep several 
numbers in mind. The variance of the noise field <r^) ~ 0.023 
due to the intrinsic ellipticities of the source galaxies is the 
same for each redshift bin; the variance of the true smoothed 
convergence field er^ is ~ 0.009, 0.017 and 0.025 for z s = 0.6, 
1 . 1 and 1 .9 respectively. The lowest threshold 2<j(N) is chosen 
so that the points considered are all in the tail of the underly- 
ing PDF: K(s)(0) > (J(s). The reason for restricting ourselves to 
the tail instead of using the full PDF is that the points which 
lie within l(T(s) of the mean probe the linear regime, where 
the shear field can be well approximated as Gaussian. In this 
regime the power spectrum, which has already been studied 
extensively, should contain essentially all of the statistical in- 
formation. 

Figure[3]shows how the fractional area statistic depends on 
the cosmological parameters f^E, 08 and wq. The data points 
in each panel are the fractional areas above different thresh- 
olds for each redshift bin in the fiducial cosmology. The lines 
are calculated by varying the parameter in question by ±20% 
(and leaving all other parameters unchanged). For reference, 
we have also shown, as the error bars, the variance measured 
directly from the simulations. The dependencies revealed on 
this figure make intuitive sense: decreasing Hue (or equiv- 
alently, increasing f2,„), increasing as, or making wq more 
negative all affect the growth of structures in the direction of 
having larger amplitudes, which increase the shear. 

Table [TJ summarizes our main results for an LSST-like sur- 
vey alone, as well as its combination with Planck. We follow 
the "self-calibration" approach, in which we assume <7(n) and 
Kbias are unknown, and fit separately in each redshift bin. To 
account for the uncertainties of the six extra nuisance parame- 
ters, instead of three S/Af bins, we utilize seven bins extending 
from v = 2.0 to v = 5.0 in increments of Av = 0.5. We consider 
two scenarios for the convergence systematics: the top half of 
Table [TJ shows a rather pessimistic projection of the survey, 
where we adopt 1% priors 14 on a(N) and Kbias, corresponding 
to the additive a nd multiplicativ e errors already achieved by 
current surveys dFuet all 12008): the bottom half of Table Q] 
represents a more optimistic projection of an LSST-like sur- 
vey, where we adopt a 0.01% prior on a^) and a 0.05% prior 
on Kbias, which is the accuracy goal of future surveys. Note 
that in the last column, we adopt a conservative approach and 
only include independent priors on H m h 2 , H^h 2 , and n s (the 
diagonal approximation instead of the full CMB Fisher ma- 
trix). These three parameters are expected to be constrained 
by Pl anck alone to 0.9%, 0.6% and 0.4% respectively dHul 
l2002h . 

As noted above, the cosmological dependence of the frac- 
tional areas, similar to galaxy cluster abundance, comes in 
through two quantities: (1) the variance of the reduced con- 
vergence, (k 2 ) /t^in' tnus most sensitive to the amplitude of 
the density fluctuation a%; (2) the minimum value of the con- 
vergence, K m ; n , thus most sensitive to the matter content fi m , 

13 www.rssd.esa.int/index. php?project=PLANCK 

14 We give equal prior to each of the redshift bin. 
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TABLE 1 

Calibrated Cosmological Parameter Constraints from LSST using the fractional area statistic and adding Planck priors: 
AS7„,/r = 0.0012, AQ b h 2 = 0.00014, and An s = 0.035. For the pessimistic scenario, we adopt 1% priors on the additive and 

MULTIPLICATIVE ERRORS ALREADY ACHIEVED BY CURRENT SURVEYS. FOR THE OPTIMISTIC SCENARIO, WE ADOPT 0.01% AND 0.05% PRIORS ON THE 
ADDITIVE AND MULTIPLICATIVE ERROR RESPECTIVELY, WHICH IS THE GOAL OF FUTURE SURVEYS. 



Parameter Constraints 


LSST (F) 


LSST {F) + Planck (priors) 


Pessimistic 






/\w., 

i_l YVp 


0.094 


0.022 


**p 


0.50 


0.94 


Awo 


0.55 


0.12 


Aw a 


1.6 


0.25 


A^DE 


0.046 


0.0095 


Acrg 


0.047 


0.0080 


Optimistic 






Aw p 


0.028 


0.012 


Zp 


0.60 


0.60 


Awo 


0.16 


0.043 


A\Va 


0.42 


0.11 


AQ DE 


0.015 


0.0038 


Acrg 


0.031 


0.0029 



or f^DE as a flat universe is assumed. However, the fractional- 
area statistic also provides good constraints on the dark energy 
equation of state, for the following reason. With Planck pri- 
ors for Cl m h 2 , Slbh 2 , and n s , which are all below the 1% level, 
an LSST-like survey can determine f^E and a% very well (to 
a few percent accuracy) by using the fractional area statis- 
tic, which in turn breaks the degeneracies in the growth of 
structure between these two parameters and the dark energy 
equation of state parameters: wo and w a . 

We follow the Da rk Energy Task Force (DETF) report 
( Albr echt et al.ll2006l) and calculate the pivot scale-factor a p , 
or equivalently the pivot redshift z p , where the dark energy 
equation of state w(a) is best constrained. In the pessimistic 
case, we find Aw p = 0.022 with a p = 0.52 (z p = 0.94) for 
LSST with Planck prior, where w p = w(a p ). This is improved 
to Aw p = 0.012 with a p = 0.63 (z P = 0.60) in the optimistic 
case. The corresponding "figure of merit" 15 is (Aw p Aw fl ) _1 ~ 
180 for the pessimistic projection, which is at the interest- 
ing level [compared t o the various stage-Ill experiments; see 
lAlbrecht etail d2006l page 77)], and 760 for the optimistic 
scenario; the best-constrained pivot point is also at somewhat 
higher redshift than other proposed probes, which gives a dif- 
ferent degeneracy direction on the wo-w a plane, so we expect 
good synergy between the fractional area statistic and other 
dark energy probes. It is also different from the main degen- 
eracy direction for CMB anisotropies, the angular diameter 
distance degeneracy, which gives a p = 0.72 (z p = 0.38) for our 
fiducial cosmology. 

In Figure |4] we show the two-dimensional marginalized er- 
ror contours of various cosmological parameters. As shown 
in the bottom panel, there is still a strong degeneracy between 
^de (or equivalently f2 m , as we have assumed a flat universe) 
and erg in the self-calibration case. Consequently, this de- 
grades the constraints of wq and w„ by large factors, as shown 
in the upper panels of Figure H] We have found that, com- 

15 Note that the DETF has included and marginalized over the curvature of 
the universe in their analysis, whereas we have assumed a flat universe. But 
we expect minor degradation of this "figure of merit" once we also include 
the Planck prior on the curvature. 



paring the pessimistic scenario with the optimistic one, there 
is respectively, a factor of 2 and 2.5 degradation for the con- 
straint on w p and w a (when Planck priors are added; the degra- 
dation is more significa nt for the convergen ce-statistic only 
case). As illustrated by iHuterer et al.l (l2006t) . the prior infor- 
mation of the intrinsic ellipticity noise can help to restore the 
dark energy constraints. We also note that the prospects for 
self-calibration of systematic errors of an LSST-like WL sur- 
vey would be further improved if the extra information from 
the tomographic power spectrum or bispectrum were added in 
the analysis. 

The numbers quoted above should be interpreted as the 
lower bound of the cosmological sensitivities one can get 
from an LSST-like survey using the statistics described (with 
Planck priors). In reality, various systematics have to be con- 
sidered, which we discuss next. 

For an LSST-like survey, most of the source galaxies will 
only have photometric redshift information. To check the 
sensitivity of the fractional area statistic to the redshift un- 
certainties, we calculated the change in the source redshift 
which would cause a shift of the fractional area by more 
than its ler error. It is estimated using Az s = (d(T) /dz s )~ l x 
y / Var[7XiJ)J, where d(T)/dz s is the derivative of the mean 
fractional area with respect to source redshift. We have found 
that the systematic error (bias) of photometric redshift of the 
source galaxies needs to be below 0.002 for low redshift and 
0.01 for z > 2, w hich is within the expecta ti ons of LSST 
see a ssumption in iMa. Hu. & Hutererl (120061) : IHuterer et al.l 
2006)]. The requirement on the photometric redshift scatter 
will not be a major problem, because many source galaxies for 
each redshift bin are needed to reconstruct the convergence 
field. 

There are also several theoretical uncertainties about the use 
of the fractional area statistic. One important issue is the the- 
oretical understanding of the one-point PDF of the smoothed 
convergence field. Similar to the case of the universal halo 
mass function, the one-point PDF has not been tested exten- 
sively for different cosmo l ogies. For instance, the expression 
given by iDas & Ostrikerl (120061) . though in excellent agree- 
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FIG. 3. — Mean fractional area (F) for different S/N threshold u, different redshift z s , and its dependence on various cosmological parameters. The symbols in 
each panel show (JF) in the fiducial cosmology (with the intrinsic ellipticity noise taken into account), with the error bars measured directly from the simulations. 
Note that the error bars are correlated, and they are also enlarged by a factor of 10 to be clearly visible. The solid (dashed) curves show results when the given 
cosmological parameter is 20% higher (lower) than its fiducial value. 



ment with N-body simulations, has only been checked for one 
fiat ACDM cosmology with Q. m = 0.3. As the above discus- 
sion indicates, the theoretical prediction for the PDF has to be 
accurate at the ~ 1 percent level. It is conceivable that the 
universality might break down at this precision. However, the 
one-point PDF is such a simple statistic and its derivation adds 
almost no extra computational cost, once WL simulations are 
made. We expect our work will inspire currently on-going or 
planned large WL simulations to obtain accurately calibrated 
formulas for the PDF. 

The baryons can cool and collapse into dense regions, 
which changes the matter distribution within the virialized ob- 
jects. The effect on the two-point statistics has been shown 
to be significant (fWh ite 2004; Zhan & Knox1l2004t Uing et alJ 
I2006t iRudd. Zentner. & Kravtsovl 120071) . and will therefore 
have to be modeled in the interpretation of future WL surveys. 
This baryonic effect would also alter the PDF tail, and thus 
the fractional area statistic, es pecially when small smooth- 
ing angles are used. However, IZentner. Rudd. & Hul (2007) 
have shown that it is possible to take into account this un- 
certainty through "self-calibration", i.e., constraining simul- 
taneously the dark energy properties and the uncertain halo 
profile due to baryonic physics, e.g., halo concentration. We 
hope to explore the utility of a similar self-calibration tech- 
nique f or the one-point statistic, using the halo model ap- 
proach dKruse & Schneid er 2000) in the future. 



The mass-sheet degeneracy of weak lensing is not an un- 
certainty for the fractional area statistic. Because the full PDF 
is measured, any offset can be determined as the mean should 
be zero. Besides, one can in pri nciple use the one-point PDF 
of the filtered tangential shear ( Valage as. Munshi. & B arber 
2005) directly instead of the convergence PDF 

It is also important to compare the fractional area statistic 
with other WL statistics. We have found that the fractional 
area statistic, in combination with CMB anisotropy measure- 
ments, may be able to reach a cosmological sensitivity closely 
approaching those from using a sample of ~ 200,000 shear- 
selected galaxy clusters (Wa ng et"aTIl2004l) . The two meth- 
ods are very similar and related in practice, as the majority 
of the hig h S/N points are candidates for being real galaxy 
clusters (iReblinskv et all Il999t I Jain & van Waerbekd 120001 
Weinberg & Kamionkowski 2003). The most important gain 
for our approach is that projection effects are not a contami- 
nation, since they depend on cosmology, and therefore serve 
as a source of information instead. The fractional area statis- 
tic does not suffer from the missing or false cluster problem, 
which is indeed the main motivation of this work. The trade- 
off, however, is that counting the number of clusters as a 
function of redshift plays a crucial role to constrain the evo- 
lution of dark energy, while the convergence field is a two- 
dimensional projection of all the structures along the line of 
sight. The tomography of the source galaxies does help, but 



10 



Wang, Haiman, and May 




0.94 



0.92 



6 0.90 



0.86 




0.66 



0.68 



0.70 



0.72 



0.74 



0.76 



0.78 



DE 



FIG. 4. — Constraints on various cosmological parameters from an LSST-like weak lensing survey. The 68% C.L. contours in all panels show reasonable 
optimistic and pessimistic estimates including systematic errors for an LSST-like survey, with or without priors from a Planck-like CMB survey. The pessimistic 
scenario corresponds to the additive and multiplicative errors already achieved by current surveys. The optimistic scenario corresponds to the accuracy goal of 
future surveys, such as LSST Different contours are for LSST alone (dotted curve) and LSST with priors from Planck (dashed curve) in the optimistic scenario; 
LSST alone (dot-dashed curve) and with Planck priors (solid curve) in the pessimistic scenario. The contours are marginalized over all other relevant parameters. 



the broad lensing kernel entangles the information from a 
wide range of redshifts; tomography is also limited by the 
number of faint galaxies detected within each redshift bin (in 
order to reconstruct a convergence field that is not shot-noise 

dominated). 

The WL p ower spectrum tomography (Song & Knox l2004t 
Zhanl 12006) and the bispectrum tomography dTakada & Jainl 
2004) have been shown to be very powerful cosmological 
probes an d sen s itive t o the evolution of dark energy. For 
example, Zhan d2006l) find that an LSST-like survey can 
deliver constraints of Awo ~ 0.16 and Aw a ~ 0.36 from 
lensing alone 16 . The fractional area statistic must include 
additional information from higher order statistics, because 
the one-point PDF is dete rmined by all orders of cumulants 
(Balian & Schaeffer 1989). Furthermore, in analogy with 
cluster counts, the fractional area statistic is helped by the ex- 
ponential sensitivity of the tail of the PDF to the cosmologi- 
cal parameters. The scale information, imprinted in different 
multipole moments of power spectrum, can be retrieved with 

16 Note these numbers are not directly comparable to ours, as the analy- 
sis in Zhan 12006) has a larger parameter set including, e.g., curvature and 
running of the spectral index; his analysis also utilizes the full Planck Fisher 
matrix instead of priors, and marginalized over 80 photometric redshift error 
parameters, which we did not do. 



different smoothing scales for the fractional area statistic: the 
linear regime of the power spectrum should be recovered by 
using various larger smoothing angular scales and the non- 
linear regime by smaller scales. However, unlike in Fourier 
space where different modes are independent, using differ- 
ent smoothing scales in real space produces mass maps that 
are correlated with each other. In this work, we only uti- 
lize one particular smoothing scale Oq = 1 arcmin, which is 
within the nonlinear regime. As mentioned earlier, in the lin- 
ear regime, most of the statistical information is already con- 
tained in the two-point correlation function. Interesting future 
work includes studying how to correctly take advantage of the 
information from the scale-dependence of the fractional area 
statistic, and its complementarity to other WL statistics. 

Finally, we note that the fractional area statistic is known as 
the first Minkowski functional. In the two-dimensional case, 
there are two other functionals: the length and the genus of 
the excursion set. These three Minkowski functionals, which 
are additive and are invariant under translations as well as 
rotations, completely characterize the morphological prop- 
erties of the high convergence regions. These functionals 
have been extensively studied for CMB as useful tests of the 
Gauss ianity of the primordial dens ity perturbation field [see, 
e.g., Winitzki & Kosowsky (1998)]. Since the convergence 
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field is definitely non-Gaussian, including the other two func- 
tional which depend on derivatives of the conv ergence field 
must yield additional co smological information (Taru va"et alj 
l2002tiGuimaraesll2002l) . 

4. CONCLUSIONS 

We have shown that, in future wide field weak gravitational 
lensing surveys, a simple one-point statistic - the total frac- 
tional area of high S/N points in the convergence field - is a 
promising probe of cosmology. It is sensitive to the total mat- 
ter content of the universe and the amplitude of the density 
fluctuations, which helps breaking the intrinsic degeneracies 
in the growth of structure between cosmological parameters, 
thus constraining the properties of dark energy. 

The main conclusion of this work is that the fractional area 
statistic provides constraints on cosmological parameters sim- 
ilar to the redshift distribution of galaxy cluster abundance, 
but without suffering from the projection effects. Indeed, the 



statistic is explicitly constructed to take advantage of projec- 
tions as part of the signal. 

We expect the fractional area statistic will help achieve the 
goal of "precision cosmology" and shed light on the mystery 
of dark energy. 
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APPENDIX 

In the appendices below, we give the expressions of the two-poin t auto-correlation function and cross-correlation function of 
the convergence fields using Limber's approximation dKaiserlll9 92). We also present a way to calculate the covariance matrix 
of the fractional area of the excursion sets, assuming the true convergence field with Gaussian smoothing can be approximated 
as a log-normal random field. The calculations below are carried out in the two-dimensional case and all the random fields are 
assumed statistically isotropic. 

A. LIMBER'S APPROXIMATION 

In this section, the two-point cross-correlation function (CCF) and auto-correlation function (ACF) of the true convergence 
field are calculated by using Limber's approximation. Note that the calculations below do not require the log-normal assumption. 

The cross-power spectrum of the convergence field on two different source planes at z, and Zj, both with Gaussia n smoothing, i s 
related to the three-dimensional matter power spectrum Ps(k) by the Fourier space analogue of Limber's equation (lKaiserlll992h : 

P(S)( k -L) = / d X 2 Ps [ — ' x ) w c( k ±9o), (Al) 

JO X \ A / 

where Xi = Xfc); W\ en (x>Xi) is given by Eq. ©; Wc(x) = exp(-x 2 /4) is the Fourier transform of the Gaussian smoothing window 
function; k± denotes the Fourier modes perpendicular to the line of sight, which are the only contribution to the projected power 
under Limber's approximation. 
The CCF is then the Fourier transform of the cross-power spectrum: 



Qs)(9) = (K(S)(8uZi)K(S)(9 2 ,Zj)) 



dX Wlen(X, XiWlen(X,Xj) 
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(A2) 



where we have changed variable k' ± = k±/x- Here A 2 (k,z) = 4irk 3 Ps(k,z)/(2ir) 3 is the dimensionless nonlinear power spectrum, 

which is calculated as in Smi th et all (12003!) : Jo denotes the Bessel function of order 0; 9 = \9\ — 62 | is the angular separation of 
two field points. 

The auto-power spectrum and the ACF of K(S) are special cases of the expressions above when i = j: 



P( S )(k±)= / d x 



(X,Xs) , 



x- 



k± 
X 



X ) W%(k ± 6 G ), 



Qs)(0)-- 



dxW len (X,Xs) 



dk'. 



A 2 (k' ±1 z) 



IJ2 



Jo(k' ± xO)mk' ± x9c) 



(A3) 



B. LOG-NORMAL RANDOM FIELD 

There is one specific class of non-Gaussian ran dom fields of theoretical interest, which can be re l ated in so me functional way to 
one or several Gaussian random fields [see, e.g.. lColes & Barrowl dl987l) : also see lBardeen et alj d 1986b and lBond~& Efstathiou 
(fl987h for a general introduction of Gaussian random field]. The useful feature of these non-Gaussian fields is that when the 
functional form of this mapping is known, the statistical properties are analytically determined from the underlying Gaussian 
fields. 

A log-normal random field is the one for which the mapping from a Gaussian random field and the inverse are given by 



K = F LN (a) = \k„ 



C/JV 



exp a w a- 



ln(l + - 



J LN 

2 

v 2 



(Bl) 
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where a is a Gaussian field with zero mean and unit variance. The minimum of the distribution K min and the variance ajj^ 
] n( I + (k 2 ) /K^ in ) are needed in order to specify k. Some relevant statistical properties of the log-normal field k are listed below: 



1 . The one-point PDF of a log-normal field k is obtained from the corresponding Gaussian PDF of a by simply changing 
variable: 

V 27r<T LN t 2(J LN J l K min|+K 



2. The two-point CCF of two log-normal fields n and k' is obtained by using the joint two-point PDF of the underlying 
Gaussian fields a and a': 

e -Qi/2 

P2[a,Oi ' \Ca(6)]dada' ' = — dada', 

Qi = l—[a 2 -2C c (9)aa' + a' 2 ], (B3) 

l-C 2 c (6) 

where Cc(9) = (a(8i)a' '(#2)) is the CCF of a and a'; 6* is the angular separation. Carrying out the average gives the relation 
between these two CCF's: 

C{6) = («(£ )«'(&)) = cxp[a LN a' LN C c (9)] - 1 , 
C c (9) = —^\n[\+C{6)l (B4) 

(T LN (J LN 

3. The relation between the two-point ACF of a log-normal field and the underlying Gaussian field is a special case of the 
expression above when k' is just k: 

C(9) = exploded M- 1, 
C G (0)=-J-ln[l+C(0)], (B5) 

where C and Cc are the two-point ACF's of k and a respectively. When 9 = 0, C(0) = (k 2 ) / K 2 in thus Cc(0) = (a 2 ) = 1 as 
expected. 

C. STATISTICAL PROPERTIES OF THE FRACTIONAL AREA 

In this section, we want to calculate the covariance matrix in Eq. (TTTb by assuming that the convergence field with Gaussian 
smoothing can be approximately described by a log-normal random field K(S)(a) = Fuj{a). To specify the function Fln, the 
parameter K m ; n is given by Eq. (O and the other parameter can be calculated using Eq. ( IA3t : 



\S)I 



■■C ( s)(0)-- 



JO J0 K 



(CI) 



as 7o(0) = 1. Note that they both depend on the source redshift z s - The CCF of the underlying Gaussian field (similarly for ACF) 
is obtained by using Eq. (IB4b . where C{9) is now C(S)(#) as given in Eq. (1A21 >. 
The mean fractional area of the excursion set E v under the log-normal assumption is given by 

(H",z s )) = -erfc (^j , (C2) 

where erfc is the complementary error function; x v = F[^(va(N)) is the effective threshold for a. 
The covariance between the fractional area with threshold of // at redshift z, and the one with threshold of v at Zj is given by 

Cov[F(»,zd,F(v,Zj)l = ((A^) - (A„) (A,)) , (C3) 

■^tot 

where the statistical average is taken by using the joint two-point PDF as given in Eq. dB3b . Because the two-point PDF only 
depends on the angular separation 9 between two different lines of sight, one spatial integral can be eliminated and gives 

Cov[F(n, zd, F(v,Zj)]= — / 2n9d9 / dn H s)(zi) / dK2(S)(Zj)P2[K\(S)(Zi),K2(S)(Zj);6]- (J r (n,Zi))(J : (i',Zj)) 

= ~^~ 2n0d9 \j x dai [ da j p ^ a i^fCG(9)]-(H^Zi))(n^,Zj))\, (C4) 



High Convergence Regions 
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where we have moved the (T) (T) term inside the integral over 9, since it has no dependence on 9; x^ and x v are the effective 
thresholds of a, and a.j respectively as defined in Eq. ( 1C21 >. 
Let us denote 

/■OO /-00 

(C5) 



{T 2 [^,Zi;v,Zj\C G (9)])= I da.il da j P 2 [a.i,a j ;C G {ff)\. 

JXy J Xv 

The special case of this expression with /x = v and i = j, had been extensively studied dKaiser|[l984t iPolitzer & Wisdfl984l : 
iJensen & Szalay 1986; Kashli nskyl 199 lh to approximate the correlation function of "biased regions", e.g., ga laxy clusters. Eval- 
uation of this integral hinges on the equality (Hamilton, G ott. & We inberg 1986; Winitzki & Kosowskyl ll998h 



d(T 2 {C c )) 
dCa 



= P 2 [x, l ,x„;C c ], 



(C6) 



which can be proven by Fourier transforming the integrand P%, then integrating over a's and Fourier transforming back. This 
first-order ordinary differential equation has a boundary condition that, when Cq = 0, the two-point PDF reduces to the product 
of two independent one-point PDF. Then the double integral can be separated and gives 

(ftfaziWfiQ)) = {H^zd){T{y,Zj)). (C7) 
Thus the term in the curly braces of Eq. (IC4b is simply 



(T2[^,KZi,Zj\C c (0)])- {Tii^^'rZi^f,®)) = I dC'P 2 (x fJ ,,x 1/ ;C') 

Jo 



(C8) 



Substituting the above expression into Eq. ( 1C41 i. and changing the order of integration by assuming Cc(9) is monotonic 17 , we 



get: 



Hot J() 



cm e\C) 

dC — pxp 



2VT Z C 2 



2(1 -C' 2 ) 



(C9) 



Here 9(C) is the inverse of the CCF: Cc(9) between a's as given in Eq. (1B4I >. The covariance reduces to the variance when /i = v 
and Zi = zf. 



1 



Hot Jo 



dC — pxp 



1+C 



(CIO) 



where 9(C) is the inverse of the ACF: Cc(9) as given in Eq. 
faster than 9~ 2 at large distances. 

We have also assumed Cc(9 m - dx ) ~ 0, i.e., the size of the field # max 
convergence field [see Eq. (ID3H . The correction term of Eq. (|C9l l is 



These integrals are finite as long as the CCF/ ACF falls off 
is large enough compared with the correlation length of the 



A(Cov) : 



Ca(9 - , dC ,l-9\C)/9. 



2 

max 



2W1-C 2 



exp 



x^ 2C x^Xv + x v 
2(1-C' 2 ) 



(Cll) 



D. RANDOM NOISE FIELD 

In the simplest case where the noise due to intrinsic ellipticity and signal are independent of each other, the one-point PDF 
of the noisy convergence field K = K(S) + K^n) is just a convolution of the log-normal and Gaussian distribution. So the mean 
fractional area of the excursion set E v with K > vagr) is 



{Hv,z s )) = 



dmw) 



2Tra 2 N) 



: exp 



(^(AQ-^bias) 

2(J In) 



1 



-erfc 



V2j' 



(Dl) 



where the effective threshold x v = F^iyagq) — Kpv)) now depends on the integral variable We define F^(k) = -00 when 
k < Kmin- We also allow a nonzero mean of the noise field Kuas as m me self-calibration case. 

Similarly for the covariance, the integral in Eq. ( IC5b will be taken over two more variables due to the noise field. The expression 
involves a five-dimensional integral which, unfortunately, cannot be simplified as what is do ne in the noise-free ca s e and is very 
messy. Thus it is useful to get an upper bound of the covariance for this general case as in Winitzki & Kosowsky ( 1998). Note 
that in the limit when x^ or x v goes to —00, (J-2) is simply (T). Since the two-point PDF P 2 is always non-negative, we get 
[HFz) < min[(^ r (/i,z,)), (T(v : Zj))]- By using this inequality, the integrand of Eq. ( IC4b then has no 9 dependence. Therefore, 



Covin^ZiXH^Zj)] < -^(^)<(l-(^)>), 



(D2) 



17 Here we have assumed that the correlation function stays positive and monotonically decreases to zero as the scale becomes large. However, the actual 
matter correlation function has to turn negative at a particular scale, before approaching zero. This is because if matter is clustered on small scales, then they 
have to be "anti-clustered" on large scales to conserve the total amount of mass. This contribution has an opposite sign, which makes the scale-dependence of the 
variance of the fractional area steeper than Poisson-scaling, i.e., the variance goes to zero faster than 9~ 2 . 
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where the subscript "<" (or ">") denotes the smaller (or larger) one between two (J F)'s. The correlation leng th 9 C is the scale 
beyond which the correlation can be neglected for physical reasons and is defined as dBond & Efstathiou 19871) 



a2 _ 



(K 2 



( K (S)) + ( K (iV)) 



c ((WK) 2 ) ((V« (5) ) 2 ) + ((V Kw ) 2 ) 
where the variances of the derivative fields are given by 



(D3) 



((Vk (s) ) 2 )=-V 2 C (5) (0)= / d X X 2 Wi n (x,Xs) 



poo 

7T / dk^AHk'^zW^xOc) 
Jo 



((V Kw ) 2 )=-V 2 C TO (0) = 



2(j2 (N) 



(D4) 



as V Jo(0) = -Jo(0). We again assume max 3> C here. 
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